Optimal design of calibration signals in space borne gravitational wave detectors 
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Future space borne gravitational wave detectors will require a precise definition of calibration 
signals to ensure the achievement of their design sensitivity. The careful design of the test signals 
plays a key role in the correct understanding and characterisation of these instruments. In that 
sense, methods achieving optimal experiment designs must be considered as complementary to the 
parameter estimation methods being used to determine the parameters describing the system. The 
relevance of experiment design is particularly significant for the LISA Pathfinder mission, which 
will spend most of its operation time performing experiments to characterise key technologies for 
future space borne gravitational wave observatories. Here we propose a framework to derive the 
optimal signals — in terms of minimum parameter uncertainty — to be injected to these instruments 
during its calibration phase. We compare our results with an alternative numerical algorithm which 
achieves an optimal input signal by iteratively improving an initial guess. We show agreement of 
both approaches when applied to the LISA Pathfinder case. 


I. INTRODUCTION 

LISA Pathfinder [? ] is an ESA mission with NASA 
contributions designed to test key technologies for the 
detection of gravitational waves in space, like the pro- 
pose eLIS A [? ] . The main scientific goal for the mission 
is expressed in terms of a differential acceleration noise 
between two test masses in nominally geodesic motion 
down to a level of Sa & = 3 x 1CV 14 m/s 2 /\/Hz at 3 mHz. 
The relevance of this requirement is not only its demand 
in terms of noise reduction but also the very low fre- 
quency measuring band, which introduces technological 
difficulties that can not be addressed by ground based 
gravitational wave detectors due to the so called seismic 
wall [? ]. 

The LISA Pathfinder mission is currently planned to 
have a six month operations period at the Lagrange 
point LI that will be split between the two experi- 
ments on-board: the European LISA Technology Package 
(LTP) and the American Disturbance Reduction System 
(DRS). This leads to a very short operation period of 
roughly three months for the complete characterisation 
and achievement of the scientific goal for the LTP. 

It is worth noticing that, after the demonstration of the 
technology readiness, a second — yet not less relevant- 
objective of the mission is a detailed characterisation of 


the noise contributions to the main scientific measure- 
ment. An extensive list of experiments has been put 
forward by the scientific team including experiments to 
characterise the optical metrology [? ] , the inertial sensor 
instrument [? ] , the effects of the thermal [? ] and mag- 
netic [? ] environment, and pure free-fall experiments 
that aim to measure acceleration noise in an configura- 
tion that is even more representative of eLIS A [? ]. All 
these runs need to be executed via tele-commands using 
a daily 8 hours communication window with the satellite. 
Internal constraints in pre-processing and validation of 
tele-commands will add a latency from 2 to 3 days be- 
tween the definition of a tele-command sequence and its 
execution on the spacecraft. 

The planning of experiments represents therefore a cru- 
cial part of the mission and needs to be optimised accord- 
ingly to make sure that the information obtained from 
each experiment is maximised. As part of this effort a 
MATLAB toolbox has been developed with the specific 
aim to deal with the LTP data during flight operations [? 
]. Among the different methods and capabilities of this 
tool, much attention has been paid to the improvement 
of the methods to obtain precise parameters from the ex- 
periments [????]. These have been tested with 
simulated data, taking into account the expected noise 
performance of the Pathfinder mission, in a series of mock 
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data challenges with data generating using the analysis 
software’s built-in modeling and simulation tools. Agree- 
ment between methods was also checked with data gener- 
ated from an independent spacecraft simulator developed 
by the prime industrial contractor, as was the case in the 
LISA Pathfinder operational exercises [? ]. 

These analyses focused on the parameter estimation 
strategy and the achievement of an optimal precision in 
the parameters obtained, following the heritage of previ- 
ous simulated data exercises, like for instance the LISA 
Mock Data Challenge [? ] that focused on the problem 
of astrophysical parameter estimation from LISA data. 
Unlike the problem of astrophysical data analysis, in the 
LISA Pathfinder case, the measured signal is the response 
of the LISA Pathfinder system to some injected input sig- 
nal that was specified by the telecommand Hie. In other 
words, there exists the opportunity in LISA Pathfinder 
to design the injected signals so that the measurement 
of the system parameters is optimised. The operators of 
ground-based gravitational-wave detectors have a similar 
opportunity to design signals when characterising the re- 
sponse of their instruments to various noise sources but, 
given their easy access to their instruments, not as much 
emphasis is placed on optimising signal injections. LISA 
Pathfinder thus represents a scenario where careful signal 
design would produce the most benefit. 

In the following we propose a general framework which 
allows the optimisation of the input signals applied to a 
given system. Experiment design can be described in 
most cases as an optimisation problem for a given figure 
of merit, which typically relates to an scalar of the Fisher 
information matrix [? ? ]. Although the description 
used here applies to a general case, in the current work 
we will be mostly interested in the application to the es- 
timation of the main parameters governing the combined 
dynamics of the test mass and the spacecraft in LISA 
Pathfinder. Hardware on-board the satellite imposes us 
a further limitation which is only to consider sinusoidal 
signals as input signals. 

This work is organised as follows. In section II we 
introduce the problem of experiment design and the no- 
tation used in this work. Section III describes a numeri- 
cal algorithm to optimise the signal to be injected given 
a model, and its application to a simple case. In sec- 
tion IV we introduce the LISA Pathfinder model used 
for our analysis in section V and finally present our con- 
clusions. 


II. FISHER MATRIX ANALYSIS 


applied to test the system and H(w; 0) is the matrix 
whose components, 7L;j(w;0), contain the transfer func- 
tion describing the dynamics of the system in the fre- 
quency domain with a dependence on a set of parameters 
0 = {0i, • • • , On}- n(u) describes the noise contribution 
of our instrument. 

The likelihood function is the probability to observe 
a measurement for a given set of parameters describing 
that system. Assuming that the data is Gaussian dis- 
tributed, the likelihood for our system will be 

p(o |0) = [27rS]- 1 / 2 exp[-i(d-H(0)-s)S” 1 (d-H(0)-s)] 

( 2 ) 

where S is the noise covariance matrix. Experiment 
design is based on the analysis of the Fisher matrix, 
whose elements are defined as 


F — 
r n 


/ / dlog(p(dj0)) \ T / <91og(p(o|fl)) 
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( 3 ) 


which can be used to set limit for expected covari- 
ance matrix of the parameters, know as the Cramer-Rao 
bound [? ] 

cov[0 i; 0j\ > F 1 (4) 


The decomposition of the Fisher matrix into eigenval- 
ues and eigenvectors will prove to be very useful in the 
following sections. Given a NxN Fisher matrix F, de- 
fined by a set of N parameters, the eigenvectors u and 
eigenvalues, A, always fulfil 


Fa = Am 


( 5 ) 


The eigenvectors can be used to diagonalize the Fisher 
matrix according to the following property 

F = R T AR (6) 

where the columns of the matrix R are the (normalised) 
eigenvectors of F and A is a diagonal matrix with the 
eigenvalues in the diagonal. Notice that R can be under- 
stood as a rotation matrix that can be used to express the 
vector of our initial parameters, 0, in the new diagonal 
basis ft, 


C = R0 (7) 

from where we obtain our new set of parameters in the 
diagonal basis, (. 


A. Definitions and notation 

In the following we will describe a given system as 

d(w) = H(w; 0) s(w) + n(w) (1) 

where o is a vector with the measurements being con- 
sidered, s is a vector with injection signals that can be 


B. Fisher matrix tomography 

To compute the Fisher matrix we need to follow Eq.(3). 
We notice though that even for this simplified problem 
the straightforward application of this expression leads 
to long expression that make difficult a further analyti- 
cal treatment. To avoid cumbersome expression as much 
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as possible we expand the Fisher matrix in its different 
composing terms. In the particular case of an experiment 
with M inputs and N outputs, we may write the elements 
of our Fisher matrix as: 

M N 

Fjj ^ ^ ) Frnnpq,ij ($) 

n,q— 1 m,p— 1 

where 

-^mnpq^ij = }mp \P@j ^PQ ( ® $q 

' . . . ' " ( 9 ) 

The definition of the Fisher matrix allows us to com- 
bine the information of different experiments by adding 
their Fisher matrices. However, in this case, we use this 
same property in the opposite direction: to split a single 
experiment as the combination of simpler independent 
experiments. This tomography from the greek tomos, 
which means “part” — will be particularly useful to inter- 
pret the Fisher matrix since we can split each experiment 
into the contribution of each transfer function and study 
them independently. The F mnpqt ij term can be under- 
stood as the mp— component of a Fisher matrix corre- 
sponding to an experiment which only considers a sinu- 
soidal input applied to the nq— channels. We notice here 
that if the noise covariance matrix, {S _1 } mp , would be 
diagonal we could consider each F mnpn ^j as the contribu- 
tion corresponding to a given transfer function H mn ((-)). 
However, cross-couplings between our channels imply a 
mixing of the different transfer function contributions. 

III. NUMERICAL DESIGN OF INPUT SIGNALS 

The analytical solution has a limited application and 
becomes unfeasible for complex systems. A more flexible 
and general approach is to iteratively optimise the input 
spectrum based on a scalar criterion. As we discussed 
above, a useful choice for this purpose is to use the de- 
terminant of the covariance matrix which is a measure 
of the volume of the uncertainty ellipsoid. For computa- 
tional simplicity, the inverse of the Fisher matrix is used 
as an approximation of the covariance matrix. Since we 
are working in a high SNR regime, it is also a good ap- 
proximation. 

The experimental design problem can be stated as how 
to choose an input signal that allows the optimisation of 
the function det[F _1 ] given the constraints of our partic- 
ular experiment. The literature suggests several options 
for a scalar figure of merit to use as a minimisation crite- 
ria including the minimisation of the trace of the inverse 
of the Fisher matrix and the maximisation of the deter- 
minant of the Fisher matrix. In our case, we will stick 
to the latter as it is the equivalent to minimising the 
uncertainty ellipsoid in the parameter space. 

This can be done analytically for simple models al- 
though no closed form exist for most problems. The 
usual strategy is to describe the problem as a numeri- 
cal minimisation problem as we show below. 


A. Dispersion function 

The dispersion function v{u) is a useful tool for exper- 
iment design. As shown below, it can be used iteratively 
to improve the input spectrum. For a given an input 
spectrum 

F 

xM = (|X(wi)| 2 . . . |X(w F )| 2 ) with \X(cj k )\ 2 = 1 

k = 1 

(10) 

the dispersion function is defined as [? ] 

v(uj) = trace [F _1 (x)F(u;)] (11) 

where F(y) is the information matrix from the power 
spectrum y(w) and F(w) is the information matrix from 
a single frequency input with normalised power spectrum 
|X(w)| 2 = 1. The dispersion function defined in this way 
has some interesting properties. In particular, it can be 
shown that the dispersion function is proportional to the 
variance of the transfer function of the system and that it 
is normalised to the number of unknown parameters, Ng. 

Given the previous properties of the dispersion func- 
tion, we can build an algorithm that minimises the co- 
variance matrix in an iterative way. The steps to do so 
are the following: 

1. Select a set of frequencies {oj\ ,...ujr} within the 
frequency band of interest and distribute the power 
equally over these frequencies. This constitutes the 
initial design. 

2. Compute the dispersion function for the F frequen- 
cies. 

3. Create a new design according to: 

X,: + i(wfc) = Xi(wfc) X Vi(u] k )/ N e (12) 

4. If rna x(is(xi, 0J k )) ~ Ng) < e for a sufficient small e, 
then the optimum design is found. If not, we return 
to step 2. 

It can be shown that the algorithm converges to an 
optimal design. [? ] 

In order to prove the efficiency of the previous numer- 
ical design method we test it in the case of a harmonic 
oscillator. We can analytically compute the Fisher ma- 
trix for this problem to obtain an expression which, as 
expected, shows a maximum of the spectrum at the natu- 
ral frequency of the oscillator, wo- This value is therefore 
the one that minimises the volume of the error ellipsoid 
in the parameter space and hence, the one that the nu- 
merical method described in the previous section should 
retrieve. 

In order to check the validity of our methodology, we 
generated a time series of 10 000 seconds of white noise 
with variance cr = 10~ 5 that we consider as our initial 
input design. We choose white noise in order to weight all 
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FIG. 1. Evolution of the algorithm to optimise the input sig- 
nal for the harmonic oscillator case. The algorithm promotes 
the natural frequency of the oscillator to o = 0.07. 


in (14) can be expressed as 2-dimensional vectors with 
components referring to the xi and xa channels respec- 
tively, 



where subindicies 1 and 2 refer to the first and sec- 
ond test mass, subindex A refers to differences between 
the first and second test mass, and capitalised subindices 
(such as force noise on the spacecraft, g n) refer to the 
spacecraft. The last equation in (14) shows how g n 
is only measured in the first channel. On the other 
hand, the differential channel is sensitive to any differ- 
ential force noises applied to the first and the second test 
masses. 


frequencies equally. We consider an harmonic oscillator 
with damping ratio £ = 0.01 and natural frequency ujq = 
0.07, and then we run the algorithm as described above. 
The result is shown in Fig. 1 where we show the evolution 
of the input signal as proposed by the algorithm. As 
shown, two iterations are enough for the algorithm to 
promote the natural frequency of the oscillator ojq among 
the others. 


IV. LISA PATHFINDER MODEL 


The matrices describing the dynamics of the LISA 
Pathfinder system are: 


D 

C 

S 


e 2 , ,2 1 mi 2 1 m 2 . ,2 
1 m sc i ra S c 1 

m 2 , ,2 
UJn 

ra sc 2 


- Wf 

s 2 + 

-Hdf 

0 ^ 


0 

H sus J ’ 


Sn 

S12 \ 


S21 

S22 ) ’ 



(14) 


In order to apply this methodology to LISA Pathfinder 
we will need first to define a model for the experiment. 
In the following we introduce the notation to describe 
the combined dynamics of the two test masses and the 
satellite required for the analysis. The same description 
with small variations can also be found in [? ? ? ]. 


A. Equation of motion 

The measurement on-board the satellite is usually ex- 
pressed as 

o = (D-S- 1 + C)- 1 (-Co I +s„+D-S- 1 o„) (13) 

where D is the dynamical matrix, C is the controller, 
and S stands for the sensing matrix, which translates 
the physical position of the test masses into the inter- 
ferometer readout, o. Subindex n stands for noise quan- 
tities, either sensing noise ( o n ) or force noise ( g n ), and 
subindex i stands for the injected signals (o»). Restrict- 
ing ourselves to linear motion along the axis between the 
two test masses (the degree of freedom that is measured 
by the interferometer), each of the dynamical variables 


where uq and io 2 are the stiffnesses — the steady force 
gradient across the test mass housing per unit mass [? ] 
— coupling the motion of each test mass to the motion 
of the spacecraft; H^f and H sus are the drag- free and sus- 
pension loops controllers, respectively. For the remainder 
of this work, it is assumed that i/df and H sus are known. 

For our current analysis we will assume some approx- 
imations in these expressions in order to keep the main 
scientific information and, at the same time, keep the 
expressions as simple as possible. For that reason, in 
the following we will eliminate the back reactions terms, 
mi = m 2 << rnsC) consider that the sensing matrix 
cross-couplings are zero Si 2 = S 21 = 0. For convenience, 
we will take the calibrations Sn = S 22 = 1. Taking into 
account these assumptions we can derive expressions for 
the transfer functions describing the system. We consider 
input signals injected at the guidance input port which 
we expressed as in Eq.(14), hence the transfer function 
is defined by 


H = (D • S -1 + C) -1 (— C) • Oi 
_ f Hn(Q) H 12 (Q) \ f On \ 
V H 2 i(Q) H 22 (Q) j ^ OiA J 


(15) 
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TABLE I. LPF noise model parameters . p 2 and pi param- 
eters correspond to frequencies in [Hz] and pi to amplitude 
spectral densities in [m/vHz] and [N/vHz] for read-out noise 
and force noise, respectively. 


Noise Parameters 

Parameter 

Onl / O n A 

9nl/ 9n2 

Qn 

Pi 

3.6 x 10~ 12 

7 x 10 -15 

2.5 x 10~ 10 

P2 

10 x 10~ 3 

5 x 10 -3 

12 x 10 -3 

P3 

4.2 

3 

3.8 

Pi 

1.8 x 10“ 3 

4 x 10" 4 

1 x 10~ 3 

Ps 

8 

8 

8 


where the transfer functions are given by 

( 16 ) 

( 17 ) 

( 18 ) 

( 19 ) 

where we realise that H12 is zero because this is propor- 
tional to the parameter S\ 2 , which we consider to be zero. 
At the same time, we see that the cross-coupling from 
drag- free to differential channel, H 2 \ 1 is proportional to 
the differential stiffness, — uj 2 . 


H u = 


ffdf 

w 2 - to 2 + H d f 


H 12 = 0 
H21 = 


H d f (w| - w i) 


H 22 = 


(w 2 - uf + H d f)(uj 2 - + H l[s ) 

Hif s 


uj 2 — uj 2 + H\f s 


B. Noise model 


Our study of the injection scheme in LISA Pathfinder 
relies on the Fisher matrix which, in turn, depends on 
the noise model used for those noise sources identified 
in Eq.( 14 ). These are: interferometer read-out noise for 
both channels — o nl and o n a — force noise applied to the 
test masses - g n 1 and g n2 — and force noise applied to the 
spacecraft — g^. Following [? ], we will characterise each 
of these with the five parameters, p\,,, 5 in the expression 


S{w) = Pi 



1 


( 


UJ 

2 TT p 2 


)P 3 


1 


( “ 

V 27 r P 4 


)P 5 


1/2 


(20) 


Applying the parameters in Table I we obtain the mod- 
els in Fig. 2 for the noise spectra of the two main interfer- 
ometer channels. We can compare the predictions from 
this simplified model to simulations coming from a de- 
tailed state-space simulator containing a much elaborate 
model of the instrument, for instance delays, actuators, 
and component noise models [? ]. As seen in Fig. 2 , 
our simple parametric model agrees well with the noise 
obtained from the state-space model. 



Frequency [Hz] 


FIG. 2. Comparison of the noise spectra for the two main in- 
terferometers channels for an analytical simple model (dashed 
line) and a noise data stream generated via a LPF state-space 
model (solid line). 


V. CALIBRATION SIGNALS FOR LISA 
PATHFINDER DYNAMICS 

During operations, LISA Pathfinder will run an ex- 
haustive characterisation campaign with the objective of 
calibrating the instrument and identifying the main noise 
contributions. Here we consider one set of experiments 
targeting the calibration of the dynamical parameters 
governing the combined motion of the two test masses 
and the satellite. For these particular set of experiments, 
the calibration procedure consists of the injection of a se- 
quence of sinusoids — the only available waveform in the 
flight software — at different frequencies at a number of 
input ports. For this work, we will focus on injection in 
one of the two main interferometer channels. However, 
the methodology can be easily applied to the remaining 
degrees of freedom. 


A. The F2121 term 

In order to demonstrate our method we consider an 
injection applied to the drag-free channel. In our frame- 
work this experiment would be completely described by 
the sum 

Fij = E (21) 

m,p 

where the indices i and j run over the parameters. The 
most general case (7 degrees of freedom) correspond to 49 
terms. This is not approachable analytically so we focus 
our attention on one term with particular relevance, the 
F2121, which can be expressed as: 

F2121,ij = {S l } 22 x 

[d e ,H 21 (u J )} T [d e] H 21 (oj)\ KM | 2 ( 22 ) 
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FIG. 3. Left: evaluation of the eigenvalue A 2 of the F 2121 term for a single frequency injection. Right: Output of the numerical 
algorithm for the optimisation of input signals based on the dispersion function applied to a LISA Pathfinder state-space model 
after 25 iterations. In this case, white noise was injected into the drag-free channel. In both cases (analytical and numerical) 
the analysis is repeated by rescaling a factor 25 and 50 the value of the second test mass stiffness, which is originally considered 
to be u>o = —22 x 10 -6 s -2 . 


This term quantifies the effect of the injection in the first 
channel as measured by the highly-sensitive differential 
channel. Under the assumptions discussed in sec. IV A, 
the only parameters that impact this term are the two 
test mass stiffnesses, which enter through the term in 
Eq. (18). Due to this simplification, we can describe this 
problem in analytical terms. Eq.(22) turns into a 2 x 
2 matrix that we can easily decompose in the related 
eigenvectors, 



' a(wi)/3(wi) ‘ 


a(w 2 )/3(w 2 ) ' 

Ui = 

a(uj 2 ) (3 (u) 2 ) 

u 2 = 

a(w i)/3(wi) 


1 


1 


1. Single tone input: undetermined solution 

If we assume that the content of our input signal is a 
sinusoid with a fixed frequency, we know from the previ- 
ous eigendecomposition that we will not be able to solve 
the problem since we have only one valid eigenvector for a 
2-dimensional problem. Nonetheless, we explore the sin- 
gle frequency solution in order to determine how much 
information can we get from the system in such a case. 

We proceed to diagonalise the Fisher matrix as in 
Eq. (6) 

F = R t AR 


where 


a(x) = H d {(uj) — x 2 + w 2 (24) 

/3(x) = H\f s (w) — x 2 + w 2 (25) 

and associated eigenvalues, 


Ai = 0 


X r^x-n n 2 , .X a 2 (wi)/3 2 (wi)a 2 (w 2 )/3 2 (a;2) 
Ai-{£ I22 HM a 4 (26) 


Since eigenvalues are directly related to the inverse of 
the expected uncertainty on the associated parameter, 
we conclude that this measurement can only constrain 
parameters in the u 2 direction while the direction u\ has 
an associated uncertainty that tends to infinity. 

It is important to notice at this point that the analysis 
we perform in the frequency domain implicitly assumes a 
unique frequency, i.e. an input signal which is a sinusoid 
at a given frequency. In the following section we explore 
which information can we obtain in such a case. 


from which we obtain a diagonal system with a unique 
eigenvalue, A 2 , given by equation Eq.(26). In Fig. 4 we 
explore this expression as a function of the frequency of 
the injection. We see that the eigenvalue has a peak when 
the input is injected at a frequency around f = 1.25 mHz. 
This becomes more evident if we increase the value of , 
as shown in the figure. 

The value of f = 1.25 mHz is therefore the best fre- 
quency for a signal composed with a unique frequency 
component for the experiment under study. Indeed, by 
maximising the Fisher matrix we are reducing the er- 
ror on the parameter space. However, it must be noted 
that this is not necessarily an optimal solution since we 
are dealing with a single frequency injection scheme that 
leads to singular Fisher matrix. 

A second consideration to take into account is that 
when diagonalising our system, our parameters are ex- 
pressed in a new basis which corresponds to applying a 
rotation matrix R to the original vector of parameters 
0 = {cof,uj 2 }. In doing so, we obtain a new set of pa- 
rameters £ = R-0. For the configuration under study the 
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FIG. 4. Determinant of the F 2121 term for an injected signal with two independent frequencies. The determinant is evaluated 
for the case = — 22 x 10 -7 s -2 (left) and cj 2 = 50 x (—22 x 10 _7 )s~ 2 (right). 


combination of parameters corresponding to the non-zero 
eigenvalue is proportional to the sum of stiffness, i.e. 

£2 oc cc 2 T CC 2 (27) 

confirming that a single frequency signal is not able to 
break the degeneracy between the two parameters in our 
system. 

We are now prepared to compare the results obtained 
analytically with the prediction of the numerical algo- 
rithm based on the dispersion function — Eq. (11). To 
do so we inject a white noise data stream to the input 
channel under consideration, which for the analysis of 
the -F2121 term is the drag-free channel. In this partic- 
ular case, we consider as our initial input a white noise 
time series of 10 5 s and cr = 10 -6 m 2 . 

In the right panel of Fig. 3 we show the resulting nor- 
malised power spectrum of the input signal as retrieved 
after 25 iterations of the numerical optimisation algo- 
rithm. The algorithm promotes the same frequencies 
that maximised the eigenvalue of the Fisher matrix F2121 
as can be seen in the left hand figure. Moreover, we per- 
formed the analysis by rescaling the value as in the 
study of the eigenvalues. Here we observe again how the 
numerical algorithm selects the f = 1.25 mHz frequency 
when approaching the case where is rescaled by a fac- 
tor of 50, proving the consistency between the analytical 
and the numerical approach. 

It is worth stressing the agreement between the two 
approaches shown, given that they are not based in the 
same description of the instrument. Whilst the analyt- 
ical derivation is funded in the expressions derived here 
from Eq. (14), the numerical approach has its roots in 
the numerical computation of the dispersion function 
Eq. (11) which uses a state space representation of LISA 


Pathfinder. The difference also lies oi how the instru- 
ment noise enter in the analysis. While we analytically 
compute the term S22 in Eq.(26), the noise enters in the 
numerical analysis through the evaluation of the Fisher 
matrix in Eq. (11). In the later, the noise spectrum of the 
instrument is computed by generating time series with 
the LPF state-space model configured with no signal in- 
jections and then computing the power spectrum. 

2 . Two tones input: full-rank solution 

Here we take advantage of the analytical solution to 
go one step further and explore the case of an input sig- 
nal composed by two sinusoids. In order to combine the 
information of more than one sinusoid frequency in the 
input signal we add the Fisher matrices corresponding 
to each frequency. Our experiment will therefore be de- 
scribed by 

AT=2 

^2121,1? = ^2121,ij(Wfc) (28) 

k 

where each ^2i2i,*j(wfc) corresponds to the contribution 
of a single sinusoid injection to the final experiment’s 
Fisher matrix. 

We first explore the rank of the F2121 matrix when eval- 
uated for different combinations of these two input fre- 
quencies. Given that F2121 depends on two parameters, 
results show that most combinations of frequencies are 
able to reach the condition rank(F 2 i2i) 7^ 2. In fact, only 
when the two frequencies are equal -and we come back to 
our previous case- we will not be in a full-rank situation. 
This allows us to go one step further and explore which 
combination of frequencies are optimal, in the sense of 
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FIG. 5. Expected error on parameters for an injection in the drag-free channel considering utf and wf as the only relevant 
parameters. Black corresponds to the initial proposal of a white noise input, blue represents the expected error for the input 
signal as obtained with the proposed numerical algorithm after 25 iterations. 


maximising the Fisher matrix, i.e. minimising the el- 
lipsoid error volume in the parameters space. Figure 4 
shows the value for the determinant of the F 2121 term 
as a function of the two injection frequencies. We ex- 
plore the determinant for two different configurations of 
the experiment: the standard with = —22 x 10 -7 s -2 
and, as before, rescaling = 50 x (—22 x 10~ 7 ) s -2 . As 
expected, the determinant shows symmetry since the two 
injection frequencies in Eq.(28) can be interchanged pro- 
ducing the same output. The determinant drops to zero 
at the diagonal since, as commented above, an injection 
with two equal frequencies sinusoid does not lead to a full 
rank solution. It is interesting to see that when we set 
50 x W2, a notch appears at the frequency f = 1.25 mHz 
that we found as a maximum in the single injection case. 

In the standard configuration, the maximum of the 
F 2121 determinant appears for frequencies in the very low 
frequency regime (f < lrnHz). If, for practical reasons, 
we set one of the two injections to be fi = 0.1 mHz the 
maximum of the function displayed in Figure 4 appears 
for a second injection at f 2 = 0.3 mHz. With these two 
values we can proceed to estimate the expected errors 
on the parameters, by evaluating the Fisher matrix in 
Eq.(28). By assuming two sinusoid injections with two 


cycles each at the obtained frequencies fi =0.1 mHz and 
fi = 0.3 mHz with and amplitude of 10~ ' m, we can eval- 
uate our expression for the Fisher matrix term F 2121 , 
obtaining a 7 x 10 -3 % and 6 x 10 -3 % relative error esti- 
mate for the two stiffness parameters w 2 and respec- 
tively. It is worth reminding here that these are optimal 
errors representing the contribution of the F 2121 term of 
the Fisher matrix to the overall experiment. We con- 
sider it as a useful example to show the capability of 
the framework here proposed to disentangle the different 
contributions to the experiment. However, the precise de- 
termination of the expected error for a given parameter 
requires the evaluation of the full Fisher matrix, which is 
composed in the analytical description of 49 components 
for the drag-free injection experiment. Hence, analysis 
considering the whole system are, in most cases, more 
suited for a numerical approach. 

In order to evaluate the improvement on the estimate 
of the parameters, we run the analysis using the numer- 
ical algorithm introduced in Sec. Ill assuming an injec- 
tion in the drag-free channel and considering only the two 
stiffness w 2 and as relevant parameters. As described 
above, the algorithm evaluates the Fisher matrix at each 
step so we can trace how the expected errors for each 
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parameter improve by modifying the input signal. The 
improvement in the error, as given by the Fisher matrix, 
is shown in Fig. 5, where we compare expected error on 
the parameters at the 1st and at the 25th iteration. The 
input signals associated with these two cases corresponds 
to a white noise injection for the first iteration that turns 
into a signal focusing all the power at f = 1.25 mHz af- 
ter 25 iterations. The results show a clear improvement 
in the expected error on the parameters which decreases 
roughly by an order of magnitude. 

VI. CONCLUSIONS 

LISA Pathfinder and future space-borne gravitational 
wave detectors will require precise calibration of their 
dynamical systems in order to operate at their design 
sensitivities. Given the operational constraints for such 
missions, the design of injection signals used for calibra- 
tion is a key aspect for efficient characterisation of the 
instrument. 

We have introduced a methodology to design experi- 
ments for these instruments based on the minimisation 
of the uncertainty ellipsoid in parameter space. This 
methodology allows one to decompose the Fisher infor- 
mation matrix in its different contributions, each related 


to a unique physical coupling — or transfer function — of 
the experiment. By studying these contributions we can 
evaluate the expected error for a given spectrum of the 
injected test signal. 

We have compared this with a numerical algorithm 
capable of generating an optimal input signal by iter- 
atively improving a proposed input spectrum. The al- 
gorithm uses the dispersion function of the system to 
promote those frequencies which minimise the error on 
the parameters under study. We have applied both tech- 
niques to one example of LISA Pathfinder injection ex- 
periments, obtaining agreement in the injection signals 
obtained with both approaches. 

As an example, we have considered the contributions 
to the expected error for a given term of the Fisher matrix 
decomposition: the F 2121 , which describes the coupling 
of the X\ (the drag-free channel) and the X12 (the dif- 
ferential channel) for the case when a signal is injected 
in the former. The analysis is however general and can 
be readily extended to other experiments within LISA 
Pathfinder. 

The methodology proposed here is general and can 
be equally applied to other instruments requiring an ac- 
curate calibration in terms of parameters uncertainties, 
such as ground-based gravitational wave detectors. 



